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ABSTRACT 

Globular Clusters (GCs) and Nuclear Star Clusters (NSCs) are typically composed by several stellar popula¬ 
tions, characterized by different chemical compositions. Different populations show different ages in NSCs but 
not necessarily in GCs. The youngest populations in NSCs appear to reside in disk-like structures, as observed 
in our Galaxy and in M31. Gas infall followed by formation of second generation (SG) stars in GCs may 
similarly form disk-like structures in the clusters nuclei. Here we explore this possibility and follow the long 
term evolution of stellar disks embedded in GCs, and study their effects on the evolution of the clusters. We 
study disks with different masses by means of detailed N-body simulations and explore their morphological 
and kinematic signatures on the GC structures. We find that as a second generation disk relaxes, the old, first 
generation, stellar population flattens and becomes more radially anisotropic, making the GC structure become 
more elliptical. The second generation stellar population is characterized by a lower velocity dispersion, and 
a higher rotational velocity, compared with the primordial older population. The strength of these kinematic 
signatures depends both on the relaxation time of the system and on the fractional mass of the second gener¬ 
ation disk. We therefore conclude that SG populations formed in flattened configurations will give rise to two 
systematic trends; (1) Positive correlation between GC ellipticity and fraction of SG population (2) Positive 
correlation between GC relaxation time and ellipticity. Thereby GC ellipticities and rotation could be related 
to the formation of SG stars and their initial configuration. 

Subject headings: methods; numerical. Galaxy; globular clusters; individual (w Centauri), Galaxy; disk, 
galaxies; nuclei 


1. INTRODUCTION 

Observations over the last two decades have shown that 
globular clusters (GCs) are not monolithic objects composed 
by a single stellar population; they host multiple stellar 
populations of stars. The different populations are charac¬ 
terized by chemical inhomogeneities in light elements, like 
He, C, N, O, E, Na, Al, Mag a nd Si, and by a strong anti - 
correlation between Na and O (iGratton et al.l 12004 2oi3)- 
Eor convenience, we will call stars with an higher content 
of He, Na and Al and a lower content of C and O “second 
generation” (SG) stars. SG stars are thought to account for a 
large fraction of the cluster total m ass, typically from ~30%- 
70% of the total number of stars dD’Antona & Caloil 120081; 
Carretta et al.||2009api|; IPTncino et al .1120101 ; iBastian & Lardol 

2015h . IBastian & Lardol (1201 5h recently argued that the 
fraction of SG stars is remarkably uniform, ~ 70%, apart 
from a few exceptions, and that the earlier reported smaller 
fractions might be considered as lower limits. The origin 
of SG stars is still debated, but they are typically thought to 
arise from gas accreted into the GC core which then forms a 
new generation of stars at its central parts. The origin of the 
gas could be either internal or external, and can originate, for 
example, from material ejected by primordial, fi rst generation 
(EG) asymptotic giant br anch (AGB) stars (IVentura et alJ 
120011; [o’ Ercole et'5l 2008 h. interacting massive EG binaries 
(deMinket^y|2009^Bashan_et alJ l2013h . very massive stars 
( Denissenkov & Hartwickl20l4l or fast rot ating massive stars 
( Decressin et al.ll2007l ; iKrause et alJ|2013ll . Material ejected 
at sufficiently low velocities (lower than the escape velocity 
from the GC) can be retained in the cluster, concentrate at 
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its center and mix with the lefto ver pristine gas and then 
fragment to form new stars (see IGratton et al.l 1201 21 and 
references therein). Different sources imply distinct He 
enrichment with respect to that of Na and a range of age 
differences between EG and SG stars, up to few hundred 
Myr in the case of the AGB polluter scenario (see e.g. 
iD’Ercole et al.ll2M^ IKrause et al.ll2013h . In order to explain 
the large SG fraction the initial c luster should have been 
initially 10-100 more massive (see IGratton et al.ll2012L and 
references therein), this issue known as the “mass budget 
problem”. 

The recently prop osed early disc accretion scenario 
(IBastian et al.l 120131) . where young low mass stars ac¬ 
crete intra-cluster gas, polluted by massive interacting 
binaries through their circumstellar discs, without requiring a 
second star formation epoch, solves the mass budget problem 
but, as any of the self-enrichement scenarios, fails to simul- 
taneously explain th e He and Na enrichment in most GCs 
(IBastian et al.l 120151) . In all the self-enrichment scenarios, 
the initial distribution of the newly formed stars st rongly 
depends on the gas configuration. Recent simulations (iBekkil 
1201(11201 ll) showed that the dissipatiye accretion of AGB 
gas, that carries a small amount of rotation inherited from 
the parent EG stars, can lead to the formation of a second 
generation gaseous disk rather than a spherical distribution 
of gas. Eragmentation of such a disk would giye rise to a 
stellar disk embedded in the spherical cluster. The interaction 
between these two stellar structures and the later cluster 
relaxation could significantly differ from that of a spherical 
cluster of a single stellar population and from that of a cluster 
wit h an embedded spherical cluster of SG stars. 

In iMastrobuono-Battisti & PereTi (120131) (hereafter Paper 
I) we began to study the eyolution of such embedded SG 
disks, and focused on the test case of an u Cen-like cluster. 
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We found that the presence of the disk leaves kinematical 
and morphological signatures on the structure of the cluster. 
These features are potentially still detectable even after 12Gyr 
of evolution. The initially spherical components is slightly 
flattened while the disk stars do not completely mix with the 
FG population. The velocity dispersions of the FG and the 
SG populations still differ after 12Gyr; the FG population 
shows a mild tangential anisotropy while the SG population 
shows significant radial anisotropy. 

Here we extend our analysis of the evolution of SG disks 
to the more general case and consider SG disks with a range 
of masses. A small spread in number ratio between the first 
and second generation stars has been recently observed (see 
e.g. iBastian & Lardol[20151) . More importantly, although 
the presence of young stellar disks in NSCs seems to be ubiq¬ 
uitous, the mass of their d ifferent stellar generations is not 
known (ICarson et alJl20L§) . Thus we explore the effects of 
these disks on the cluster evolution, and the dependence of the 
cluster properties on the fractions of SG stars in the disk. We 
also compare these results to the case of SG stars formed in 
an embedded spherical sub-cluster rather than a disk configu¬ 
ration. 

We note that recently iHenault-Brunet et'di (1201 5h have 
also studied SG formation scenarios and considered forma¬ 
tion in SG disks. When relevant we compare our results with 
their findings. 

In Section |2] we describe the methods used and the sim¬ 
ulations run. The results of these simulations are shown in 
Section 0 their observable implications are and discussed in 
Section m followed by the summary and conclusions in Sec¬ 
tion |5] 

2. METHODS 

The evolution of massive second generation disks in GCs 
has been explored by means of high-precision N-body simu¬ 
lations. In the following we describe the methods and tools 
used (similar to those used in Paper I), as well as the various 
initial conditions (summarized in Table [T]). 

2.1. The N-Body Simulations 

We run several N-body simulations of stellar disks 
embedded in an i nitially spherical cluster, using 0GPU 
(iBerczik et alJ20in) . a direct N-body code running on graphic 
processing units (GPUs). The code has been developed based 
on 0GRAPE, a previous code optimized to use GRAPE6 
cards as hardware accelerators. 0GPU code is fully paral¬ 
lelized using the MPI and CUBA C, the language developed 
to program NVidia GPUs. 

We run our simulations both on the cluster Tamnun at the 
Technion and on the cluster Cytera in Cyprus. Tamnun is 
equipped with 4 nodes with one Nvidia Tesla M2090 each 
and other four nodes with three Tesla K20m. Cytera is an 
HPC cluster with 18 GPU nodes with dual NVidia M2070 
GPUs. In order to smooth close encounters between particles 
we used a softening length of 10“^ pc. As shown in Paper I 
even a softening length as big as 5 x 10“^ pc does not affect 
the results significantly. The relative energy variations at the 
end of the simulations are smaller than AEjE = 0.01. 

2.2. Initial Conditions 

The initial spherical, EG, component of the cluster is a sin¬ 
gle mass King model , similar to the best fit model of the den¬ 
sity profile of u Cen (lMevlanll987h . with Wo = 6, core radius 


TABLE 1 

Number of particles in the FG and SG populations and relative initial 

HALE-MASS RELAXATION TIMES. 


ID 

Nfg 

Nsg 

La.o (Gyr) 

S19 

50000 

11175 

7.0 

S24 

50000 

16121 

7.0 

S29 

50000 

19193 

8.2 

S33 

50000 

25069 

8.5 

S50 

50000 

45865 

8.3 

S50s 

50000 

50000 

8.6 


Kc - 4.4 pc and tidal radius r, - 80 pc. This system is sampled 
using N - 50000 particles, whose mass is 20 Mq, accounting 
for a total mass of M = 10® Mq. 

We embed each of the SG disks we consider in the initially 
spherical FG stellar population. We assume that the cluster 
has already lost a large fraction of FG stars and we focus on 
its the long term evolution. The disk is allowed to relax in¬ 
side the cluster until it reaches a quasi-stable state, which is 
then considered to be the initial configuration of the system. 
During the pre-initial stage a fraction of the SG stars in the 
disk are expelled from the system, and are removed from the 
N-body simulation, leaving behind a disk whose mass corre¬ 
sponds to 19%, 24%, 29%, 33% and 50% of the total mass of 
the system (first-i-second generation populations), in the dif¬ 
ferent models considered (see Table[T]for more details and for 
the simulations IDs used throughout the text). The profile 
of the FG component, after the initial relaxation phase, is still 
well described by the initial parameters. At the same time 
the disk inflates achieving a flattening q=0.4; 90% of its to¬ 
tal mass is inside the core radius of the FG spherical cluster. 
The initial half-mass relaxation time, f^/Q, of each simulated 
system is also given in Table[T] T he FG spheric al component 
needs to rotate to form a SG disk (lBekkill2010l) . However, in 
our case, the initial disks have a maximum rotational velocity 
of ~ 7km/s. Since this value corresponds to a peak rotational 
velocity smaller than 0.5km/s for the FG population (see Fig¬ 
ure 3 in we neglect this initial effect. N-body 

simulation of GC systems are computationally expensive, and 
we therefore use a smaller number of particles to represent 
the FG cluster and SG disk stars, as to allow for faster com¬ 
putation. The mass of each N-body particle is therefore larger 
than the typical average stellar mass in GCs. Since this choice 
affects the dynamical time-scales, we rescale the simulation 
time to the relaxation time of the corresponding real system, 
using the procedure described in Paper I, assuming an average 
stellar mass of (m.) = 0.5 Mq for the stars in a realistic GC. 
Each system is evolved for 12Gyr of rescaled simulation time, 
comparable to the age of GCs. 

The dynamics of stars born in a disk can give rise to differ¬ 
ent signatures observable in the cluster structure, compared 
with the effects of SG stars born in a spherical sub-cluster. In 
order to compare these two different initial configurations, we 
run an additional simulation where the SG stars are distributed 
inside the core radius of the FG population with a King pro¬ 
file of the same parameters as used for the FG stars. The latter 

^ trh has been evaluated using the definition given by 4SDitzeiiri987li . esti¬ 
mating the number density and the root mean square velocity. This expres¬ 
sion of the relaxation time is in principle only valid for spherical and isotropic 
systems. For a system with coh erent velocit ies (such as a disk), the relaxation 
time would be shorter (see e.g. lHaasll2014l) . The half-mass relaxation times 
listed in Table[T]may therefore be considered as upper limits for systems with 
an embedded disk. 
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Fig. 1.— The evolution of the SG simulated disks with time. The isodensity contours ai'e given at 0, 3, 6, 9 and 12 Gyr (going from the left to the right) and for 
the 19%, 24%, 29%, 33%, 50% mass fraction disks (going from top to bottom). The density contour levels are 50, 100, 200, 300, 500 and 1000 Mo/pc^. The z 
axis is parallel to the L, component of the disk angular momentum. The disks cleai'ly inflate with time, becoming almost spherical after 12Gyr. Different systems 
have different relaxation times so their final configuration slightly depends also on this factor. 


case is simulated only for a SG population comprising 50% of 
the total mass of the cluster (including both EG and SG stars). 

The clusters are simulated in isolation; tidal effects due to 
the galaxy significantly affect only the most external regions 
of the cluster, and are neglected her43- We did not take into 
account stellar evolution. Generally, mass loss would slow 
down relaxation, and stronger kinematic signatures should be 
left. We did not consider a mass spectrum for the stars, which 

^ In order to check that the effects of an external potential can be neglected 
we run an additional simulation with the cluster orbiting a point-mass galaxy 
with the same mass of the Milky Way within 8.5 kpc, the cluster being on 
a circular orbit at this radius. The final system (we ran only the case with a 
mass ratio of 33% between the two populations) shows the same final features 
of its counterpart evolved in isolation. 


could also affect the relaxation time, however, massive stars 
are short lived (and stellar black holes are rare) and most of 
the cluster evolution would be dominated by the dynamics of 
low mass stars, for which the range of masses is modest. 

3. RESULTS 

We analyzed the morphological and kinematical properties 
of the system after 12Gyr in order to find how the initial pres¬ 
ence of a disk modifies the structure of the whole cluster, leav¬ 
ing behind long term observable signatures. 

3.1. Morphological properties 

Eigure [T] shows the isodensity contours of each stellar disk 
at the beginning of the simulation and after 3, 6, 9 and 12Gyr. 
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(1201 5h . who attributed it to the initially more centrally con¬ 
centrated configuration of the SG stars. In addition, we find 
that the higher radial anisotropy is correlated with higher mass 
ratio between the SG and FG stellar populations. This cor¬ 
relation is due to the higher density of the disk and shorter 
relaxation time in systems with a more massive SG in respect 
oftheFG. 

The right panel of Figure |4] shows the azimuthal velocity 
vg of the FG and SG populations as a function of the radius. 
In all the cases the FG populations show a mild central ro¬ 
tation with velocity between 0.5 and 1 km/s. The SG popu¬ 
lations, regardless of the mass of the disk, rotate with vg be¬ 
tween 0.5 and 2km/s (within the central ~ 20pc). The rotation 
initially increases, up to 3 pc, and then decreases with the 
distance from the center of the cluster. As expected (see also 
iFlenault-Brunet et~d1l2015h . the spherical case does not show 
any significant rotation. 


Fig. 2.— Axial ratios of both first and second generation populations after 
12Gyr of evolution. 

As already found in Paper I and by IVesperini et al.l (120131) . 
the second generation stars always remain confined within the 
central 10-20pc without completely mixing with the FG pop¬ 
ulation. The 19, 24 and 29% disks, at the end of the simula¬ 
tion, are almost spherical while the 33% and the 50% disks 
are still significantly flattened. We also note that the 50% disk 
is unstable and forms a central bulge that could lead to a faster 
angular momentum exchange with the FG population through 
collective effects rather than two-body relaxation (see also 
Paper I). 

Figure |2] shows the axial ratios of the whole system at the end 
of each simulation. The 19% disk is slightly prolate, with 
all axial ratios ~ 1, while the other systems are oblate. The 
cja axial ratio is smaller for larger disks, the 50% disk is al¬ 
most perfectly oblate while the clusters with lower mass disks 
are more triaxial. The cluster is flattened at any radius, but 
the effect is larger in the central 10-20pc, where most of the 
disk stars reside. In contrast, the spherical SG system (model 
S50s) and its host clusters system do not show any significant 
flattening, as expected. 

3.2. Kinematic properties 

Figure [3] shows the final tangential and radial velocity dis¬ 
persions of all the initially spherical components and their 
comparison to embedded disk cases. Both the disk and the 
cluster (FG) stars are centrally isotropic. 

In every case the disk stars show a large radial anisotropy 
that increases with the distance from the center (see left panel 
of Figure H]). This is true also for the FG population that, 
nevertheless, shows a milder radial anisotropy. The velocity 
dispersion of the SG population remains lower than that of the 
FG population even after 12 Gyr, as also found for the low 
mass disk simulated in Paper I. This difference is more ap¬ 
parent in the intermediate/external regions of the system. The 
level of radial anisotropy of the disk stars does not show a 
clear dependence on the initial disk mass. In contrast, as also 
seen in Figure [3] the FG population shows a distinct depen¬ 
dence of its radial anisotropy on the initial disk mass. 

The S50s case is highly radially anisotropic (see Figures 
[3] and H]). Both the first and second generation populations 
are more radially anisotropic than the systems hosting an ini¬ 
tially flattened second stellar generatio n. This induced radial 
anisotropy has been similarly found by iFlenault-Brunet et"^ 


3.3. Projected axial ratios 

In order to directly compare the morphological properties 
of the simulated systems to those of their observed counter- 
pa rts we evaluated the proje cted axial ratio qp of each system 
as (iChakrabartv et al.ll200^ 


q 


2 

P 


2 2 


cos^(i) + q^sin^(i) 


( 1 ) 


where qi and q 2 are the 3D major and minor axial ratios 
and i is the inclination of the line of sight in respect to the 
z axis. We evaluated qp for all the systems after 12Gyr of 
evolution (Figure |3. We adopted four different inclinations: 
0°, 30°, 60° and 90°. We find that the projected flattening 
grows with the disk mass. Note, however, that the observed 
flattening strongly depends on the inclination angle in respect 
to the line of sight; in particular, even very flattened clusters 
could appear spherical if seen face-on. Inclinations of 60° 
and 90° (edge-on) amplify the flattening compared with the 
other cases. 

As already found for the spatial axial ratios, the projected 
axial ratios of the cluster hosting a spherically distributed SG 
population are of order of unity and very similar, regardless 
of the adopted viewing angle. 


4. OBSERVABLE IMPLICATIONS 

In the following we consider the observable signatures of 
second generation populations and their initial configurations 
on the current structures and properties of Galactic GCs. In 
particular we focus on specific trends expected as a function 
of the fraction of SG stars and the flattened vs. spherical con¬ 
figurations. 


4.1. Flattening 

Typical GCs are not spherical but are slightly ellip¬ 
soidal (see e.g. i Gever et ^ I1983I IWhite & Shawil 119871 : 
iMackev & van den BerghI 120051 : IChen & ChenI 1201 Oh . Sev¬ 
eral origins for the observed ellipticity has been suggested, 
including primordial in ternal rotatio n , exter n al tides, and/or 
press u re anisotropy dFiestas_et^lJ_|200^; Ivan den BerghI 


120081: IVarri & BertinI 1^121 : iBianchini et al. 120 1311) ^ Mere - 
ers of globular clusters, that have been proposed as a 
possible expl anation for the prese nc e of multiple stellar 
generations (Ivan den Ber^ 119961 : ICarretta et alJ 120101 : 
lAmaro-Seoane et alJ 2013l k co uld also lead to significant 
flattening (iMakino et alJ 119911) . However, the absence of 
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Fig. 3. — Radial (solid line) and tangential (dashed line) velocity dispersions for the different simulated systems, after 12Gyr of evolution. The black lines are 
for the FG population (cluster) and the blue ones are for the SG population (disk). 




Fig. 4.— The final velocity anisotropy parameter (left panel) and the azimuthal velocity (right panel) for all the systems. The solid lines ai'e for the initially 
spherical component, the dashed lines are instead for the flattened or spherical second generation population. 


binary Galactic GCs (see e.g. linnanen et al.l Il972t ISurdinI 
ll99Ulvan den Berghl2008l) suggests that this process does not 
play an importan t role in the eyolution of GCs in our Galaxy. 
lAkiyamj (1199 Ih found that the external regions become 
more spherical with time as a consequence of grayothermal 
contraction and tidal stripping. At the same time, the tidal 
field can st retch the c luster s, making them more elongate. 
Moreoyer, I Schmidt! (119961) found that moderate bright 
Galactic GCs are rounder than fainter and brighter clusters. 
The existence of a correlation between age and flattening 
(iMakino et alJ 1199 IL and reference therein) or between 
flattening and the prese nce of bluer horizontal branches in 
clusters are still debated (lNorrislll983l : lyan den Berghll2008h . 


As also shown in Paper I, in our simulated clusters, the flat¬ 
tening is caused by the angular momentum exchange between 
the second and first generation stars. When both the first and 
second generation are spherically symmetric, the system al¬ 
ways keeps its initial morphology and does not rotate (see 
Sections 13.11 and 14.3b . Note that since it does show signifi¬ 
cant radial anisotropy, the y eloci ty a niso tropies are not linked 
with the flattening (Sections l3.2l and l4r2l) . 

Since the central density and the angular momentum in¬ 
crease with SG disk mass, both the total amount of the an¬ 
gular momentum transferred and the efficiency of the angular 
momentum exchange between the FG and SG populations are 
expected to be larger for more massiye SG disks, as indeed 
obseryed in our simulations. 
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Fig. 5. — Projected axial ratios of the simulated systems at the end of the 
runs. 

The 3D flattening is always larger than the projected one 
(see Figures |2] and |5]l, thus observations can only provide a 
lower limit for this quantity. However, statistical analysis 
and future observations, e.g potentially provided by the GAIA 
mission, could give more detailed information on the spatial 
structure of GCs, allowing a quantitative comparison between 
the theoretical expectations and the observations. 

The projected axial ratios (see Equation ([T]i) of the simu¬ 
lated systems, averaged over the radial extent and evaluated 
at the half mass radius, are shown in Figure |6] as function of 
the SG mass fraction of the total mass of each system and 
for four different projection angles (0°,30°, 60°, 90°). If seen 
face-on (0°) all the clusters look almost spherical. Different 
inclinations produce different values of the flattening, and this 
could be connected to the large spread in qp for the observed 
Galactic GCs (see Figure |6]l. The trend of ellipticity correla¬ 
tion with SG mass fraction can be clearly seen. 

4.2. Velocity anisotropies: a tracer for the mass of second 
generation stellar population 

Both the cluster hosting a spherical SG population and the 
clusters whose SG stellar populations have initially flattened 
structures show radially anisotropic profiles. Therefore, this 
feature canno t be used as a tracer for the i nitial SG configura¬ 
tion (see also iHenault-Brunet et~ani2015h . however it can be 
used as an indicator of the total mass of the SG population. 
The simulated clusters have similar relaxation times and their 
radial anisotropy grows with time. After 12Gyr of evolution, 
the larger the second generation mass (independent of the ini¬ 
tial configuration), the larger is the radial anisotropy of the FG 
population. The SG population is more radially anisotropic 
than the FG one. 

One difference between the spherical SG case and the flat¬ 
tened SG populations is that even when the mass of the spher¬ 
ical SG model is similar to that of the most massive simulated 
disk, its velocity dispersion is more radially biased than that 
found for the S50 case. This is possibly caused by the faster 
relaxation in the case of an initially flattened distribution for 
the SG population. 0 

^ If a disk-like structure relaxes faster than a spherical one (see e.g. IHaa.sl 


4.3. Rotation: the disk signature 

A bulk average rotational velocit y has been already ob¬ 
serv ed in few Galact i c GCs . Recentlv lBianchini et al.l(l2013al) 
and iKacharov et alJ (1201 4ft found that the rotation plays an 
important role in the shaping the structure of co Cen, Ml5, 47 
Tuc and NGC 4372 and it gives rise to t he flat tening of these 
clusters. Furthermore, iFabricius et al.l (1201 4ft found a tight 
correlation between rotation and flattening for the 11 Milky 
Way GCs in their sample. 

The FG population in our simulation rotates with vg smaller 
than Ikm/s. The SG population have a maximum vg between 
Ikm/s and 2km/s, and all the systems reach the maximum ve¬ 
locity at ~3pc from the center of the cluster, regardless of the 
mass of the disk. The rotation is apparent even beyond lOpc 
- the radius containing most of the second generation stars. 
As a consequence of how the initial conditions are built the 
spherical S50s system does not show any significant rotation 
in any of the components. Thus, in the AGB model or in any 
model where the SG stars form in a disk, the rotation can be 
potentially explained as the result of angular momentum ex¬ 
change between an initially flattened and slowly rotating SG 
stellar population and an initially spherical FG population. 

4.4. Core collapse 

The strength of the signatures left by the SG disk depends 
both on the mass of the disk and on the relaxation time of the 
system, but all the signatures are evident even after 12Gyr in 
any of the simulated cases. All the systems studied have re¬ 
laxation times shorter than a Hubble time and the presence of 
disks induce a more rapid core collapse of the cluster; much 
faster than that expected from non-centrally concentrated SG 
population. More massive disks lead to faster core collapse, 
and we therefore expect clusters with a larger second genera¬ 
tion fraction to be more dynamically evolved than their coun¬ 
terpart similar systems with smaller second generation frac¬ 
tions. 

5. SUMMARY AND CONCLUSIONS 

We have explored the long term evolution of massive sec¬ 
ond generation disks embedded in an initially spherical clus¬ 
ter composed of primordial, first generation stars. We have 
also compared these results with the case of SG stars in a 
spherical configuration. The initial existence of a SG disk 
gives rise to kinematic and structural signatures imprinted in 
the cluster properties. These signatures are potentially ob¬ 
servable even today after 12Gyr of evolution, and the strength 
of these signatures is correlated with the SG population mass 
fraction. These signatures become weaker and less observable 
due to relaxation processes, and should therefore be more pro¬ 
nounced in clusters with longer relaxation times; nevertheless, 
observable signatures can still linger even after more than a 
relaxation time. 

Our main findings can be summarized as follows 

• Similarly to what found in Paper I and in 
IVesperini et al.l (1201 3ft . second generation stars remain 
confined in the central 10-20pc and never completely 
mix with FG population. This is still true even after 
more than one (initial) relaxation time of the system. 

|20T1) . the radial anisotropy due to the initially more concentrated configu¬ 
ration is expected to be erased earlier in time for an initially flattened SG 
population. 
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Fig. 6.— Average (left panel) and half mass radius (right panel) projected axial ratios of the simulated and observed GCs plotted against the SG mass fraction 
in respect to the total mass of the system. 


• The evolution of the disk inside the cluster leads to a 
significant flattening of the system (see Paper I as well). 
The amount of flattening increases with the mass of the 
SG disk. 

• The least massive system becomes slightly prolate 
while the other clusters are significantly oblate, espe¬ 
cially in their central regions. However, the more mas¬ 
sive systems are oblate up to their tidal radius. There¬ 
fore, future observations of the external, less crowded, 
regions of GCs may provide us with important clues 
regarding the early star formation history in these clus¬ 
ters. 

• The most massive system is also the most flattened one 
(see Figure|2]i. Its cja ratio is consistent with tha t of the 
most flattened GCs observed in the Milky Way (iHarrisI 

[HM). 

• The projected axial ratios strongly depend on the view¬ 
ing angle and increase with the mass of the SG disk. 

• The spherical component, as well as the whole system, 
bec omes radially anisotropic a s the disk inflates (see 
also lHenault-Brunet et al.l2015h . In particular, the more 
massive the disk is, the more radially anisotropic the 
cluster becomes. 


the second generations and not of their initial configu¬ 
ration. 

• Confir ming what found by iHenault-Brunet et"^ 
(1201 5h . the rotational features of the two populations 
are remarkably distinct. While the EG population does 
not rotate, at least outside the central few parsecs, the 
SG population does show rotational features, with a 
peak in the azimuthal velocity around 3pc from the 
cluster center. The more massive the SG population 
the larger is the overall cluster ellipticity. Future 
observations of differential rotation between first and 
second generation stars could then tell us about the 
initial configuration of the second stellar generation, 
shedding light on its formation mechanism. 

• The presence of a central denser population speeds up 
the core collapse. 

In conclusion, the contemporary observation of flatten¬ 
ing, lower velocity dispersion for the SG stars, radial 
anisotropy and SG rotation could point to an initially 
disk-like configuration of the younger, light-element 
enriched populations. The strength of these signatures 
could provide information on the dynamical age of the 
clusters and also on the fraction of SG stars present in 
the system. 


• As already observed in Paper I, the EG population 
shows a larger velocity dispersion compared with the 
SG population. 

• This is true also when the SG has an initially spherical 
configuration. The anisotropy is a tracer of the mass of 


Akiyama, K. 1991, Earth Moon and Planets, 54, 203 
Amaro-Seoane, R, Konstantinidis, S., Brem, R, & Catelan, M. 2013, 
MNRAS, 435, 809 

Bastian, N., Cabrera-Ziri, L, & Salai'is, M. 2015, MNRAS, 449, 3333 


The authors wish to thank the referee for helpful comments 
and suggestions. This research was supported by the T 
CORE Program of the Planning and Budgeting Committee 
and The Israel Science Foundation grant 1829/12. We ac¬ 
knowledge the Cyprus Institute under the LinkSCEEM/Cy- 
Tera project. 


Bastian, N., Lamers, H. J. G. L. M., de Mink, S. E., et al. 2013, MNRAS, 
436, 2398 

Bastian, N. & Lardo, C. 2015, MNRAS, 453, 357 
Bekki, K. 2010, ApJL, 724, L99 





















8 


Mastrobuono-Battisti & Perets 


Bekki, K. 2011, MNRAS, 412, 2241 

Berczik, R, Nitadori, K., Zhong, S., et al. 2011, in International conference 
on High Performance Computing, Kyiv, Ukraine, October 8-10, 2011., 
8-18 

Bianchini, P., Varri, A. L., Berlin, G., & Zocchi, A. 2013a, ApJ, 772, 67 
Bianchini, P, Varri, A. L., Berlin, G., & Zocchi, A. 2013b, MmSAI, 84, 183 
Carretta, E., Bragaglia, A., Gratton, R., & Lucatello, S. 2009a, A&A, 505, 
139 

Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2009b, A&A, 505, 117 
Carretta, E., Gratton, R. G., Lucatello, S., et al. 2010, ApJL, 722, LI 
Carson, D. J., Barth, A. J., Seth, A. C., et al. 2015, AJ, 149, 170 
Chakrabarty, D., de Filippis, E., & Russell, H. 2008, A&A, 487, 75 
Chen, C. W. & Chen, W. P. 2010, ApJ, 721, 1790 
D’Antona, E. & Caloi, V. 2008, MNRAS, 390, 693 

de Mink, S. E., Pols, O. R., Langer, N., & Izzard, R. G. 2009, A&A, 507, LI 
Decressin, T., Meynet, G., Charbonnel, C., Prantzos, N., & Ekstrom, S. 

2007, A&A, 464, 1029 

Denissenkov, P. A. & Hartwick, E. D. A. 2014, MNRAS, 437, L21 
D’Ercole, A., D’Antona, E, Ventura, P, Vesperini, E., & McMillan, S. L. W. 
2010, MNRAS, 407, 854 

D’Ercole, A., Vesperini, E., D’Antona, R, McMillan, S. L. W., & Recchi, S. 

2008, MNRAS, 391,825 

Fabricius, M. H., Noyola, E., Rukdee, S., et al. 2014, ApJL, 787, L26 
Fiestas, J., Spurzem, R., & Kim, E. 2006, MNRAS, 373, 677 
Geyer, E. H., Nelles, B., & Hopp, U. 1983, A&A, 125, 359 
Gratton, R., Sneden, C., & Can'etta, E. 2004, ARA&A, 42, 385 
Gratton, R. G., CaiTetta, E., & Bragaglia, A. 2012, A&ARr, 20, 50 
Haas, J. 2014, PhD thesis, Charles University in Prague, Faculty of 
Mathematics and Physics, Astronomical Institute, V Holesovickach 2, 

18000 Prague, Czech Republic 


Harris, "W. E. 1996, AJ, 112, 1487 

Henault-Brunet, V., Gieles, M., Agertz, O., & Read, J. I. 2015, MNRAS, 
450, 1164 

Innanen, K. A., Wright, A. E., House, E. C., & Keenan, D. 1972, MNRAS, 
160, 249 

Kacharov, N., Bianchini, P, Koch, A., et al. 2014, A&A, 567, A69 
Krause, M., Charbonnel, C., Decressin, T., Meynet, G., & Prantzos, N. 2013, 
A&A, 552, A121 

Mackey, A. D. & van den Bergh, S. 2005, MNRAS, 360, 631 
Makino, J., Akiyama, K., & Sugimoto, D. 1991, Ap&SS, 185, 63 
Mastrobuono-Battisti, A. & Perets, H. B. 2013, ApJ, 779, 85 
Meylan, G. 1987, A&A, 184, 144 
Norris, J. 1983, ApJ, 272, 245 

Pancino, E., Rejkuba, M., Zoccali, M., & Can-era, R. 2010, A&A, 524, A44 

Schmidt, K.-H. 1996, Astronomische Nachrichten, 317, 333 

Spitzer, L. 1987, Dynamical evolution of globular clusters 

Surdin, V. G. 1991, Ap&SS, 183, 129 

van den Bergh, S. 1996, ApJL, 471, L31 

van den Bergh, S. 2008, AJ, 135, 1731 

Van-i, A. L. & Berlin, G. 2012, A&A, 540, A94 

Ventura, P, D’Antona, R, Mazzitelli, L, & Gratton, R. 2001, ApJL, 550, L65 
Vesperini, E., McMillan, S. L. W., D’Antona, R, & D’Ercole, A. 2013, 
MNRAS, 429, 1913 

White, R. E. & Shawl, S. J. 1987, ApJ, 317, 246 


